September 1, 2008 15:30 WSPC - Proceedings Trim Size: 9in x 6in texte 



1 



00 

o 

O 
(N 

Local control of area-preserving maps 

'. 

. C. Chandre and M. Vittot 



Q 
U 
d 



Centre de Physique Theorique* , CNRS Luminy, Case 907, F-13288 Marseille Cedex 

9, France 

G. Ciraolo 

MSNM-CP"^ , IMT La Jetee, Technople de Chateau Gombert, F-13451 Marseille Cedex 

20, France. 



We present a method of control of chaos in area-preserving maps. This 
method gives an exphcit expression of a control term which is added to a 
given area-preserving map. The resulting controlled map which is a small and 
^ 1 suitable modification of the original map, is again area-preserving and has an 

' invariant curve whose equation is explicitly known. 

cn ■ 

■ 1. Introduction 

On ' Chaotic transport arises naturally in Hamiltonian systems with mixed 

. phase space. Achieving the control of these systems by restoring local con- 

I served quantities is a long standing and crucial problem in many branches 

of physics (in particular, in plasma physics and fluid dynamics). A method 
for controlling continuous Hamiltonian flows has been developed based on 
' the following idea: to find a small control term / for the perturbed Hamil- 

■ tonian H ~ Hq + V (where Hq is integrable), in order to have a more 
regular dynamics for the controlled Hamiltonian He = Hq + V + f . Two 
approaches have been developed : A global control aims at making the 
controlled Hamiltonian He integrable; A local control restores a particular 
invariant torus (local integrability) . Both approaches give a control term of 
order 
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Let us stress that this method of control differs from other methods 
by the fact that the controlled dynamics is Hamiltonian : This makes it 
relevant to the control of inherently Hamiltonian systems such as beams of 
charged test particles in electrostatic waves, two-dimensional Euler flows or 
the geometry of magnetic field lines. 

These two control methods have been developed for continuous time 
flows. ^'^ In Ref.,^ the global control for symplectic maps has been proposed. 
In this article, we explicit the local control method for symplectic maps. 

In Sec. 2, we derive the expression of the control term for area-preserving 
maps generated by a generating function in mixed coordinates. In Sec. 3, 
we apply the local control of area-preserving maps to two examples : the 
standard map and the tokamap. 

2. Derivation of the control term 

We consider two-dimensional symplectic maps {A, (p) i— > {A', ip') = F{A, ip) 
on the cylinder R x T which arc e-closc to intcigrability.^ In this section, our 
aim is to find a small control term / siich that the controlled symplectic map 
F + f has an invariant curve. We consider area-preserving maps obtained 
from a generating function of the form 

S{A', if) = A'^ + H{A') + eV{A', ip). 

The map reads 

A = A' + ed^V{A',^), 

ip' = ^ + H'{A')+edAV{A',ip). 

Here dAV{A','p) denotes the partial derivative of V with respect to the 
action (first variable) and d^V{A',{p) denotes the partial derivative of V 
with respect to the angle (second variable). 

We expand the map around a given value of the action denoted K. The 
generating function after the translation is : 

S{A', ^) = + H{K + A') + eV{K + A', ^). 

We rewrite the generating function as : 

S{A', ^p) = A'^ + u)A' + ev{<p) + w{A', ip), (1) 

where 

u = H'{K), (2) 
v{^) = V{K,ip), (3) 
w{A', ip) = H{K + A') - H{K) - ojA' + eV{K + A', ip) - eV{K, (p).{4) 
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We notice that w(0, (p) = for all G T. Without loss of generality, 
we assume that J^^ v{tp)dip = 0. Our aim is to modify the generating 
function with a control term / of order such that the controlled map has 
an invariant curve around A' = Q. We consider the controlled generating 
function 

Sc{A', If) = A'if + (jA' + sv{tp) + w{A', ip) + f{<p), 

where we notice that the control term / we construct does only depend on 
the angle if. The controlled map is given by 

A = A' + ev'{^)+dMA',^)+f{ip), (5) 
(p' ^ ip + Lj + dAw{A',tp). (6) 

We perform a change of coordinates generated by 

X{Ao,(p) ^ Aoip + exi(p), 

which maps {A,(f) into {Ao,(fo), and {A',(f') into {A'qjIp'q). The mapping 
becomes 

Ao^A'o + s ix'i^'o) - x'i^o) + v'iifo)) + d^w (A'q + ex'K), <Po) + /'(v'o), 
ip'o = ^o+t^ + dAw{A'o+sx'i(p'o),<Po)- (7) 
We choose the function x such that 

X{V> + i^) - Xif) = -vif)- 
By expanding v in Fourier series, i.e. v{(p) = J2kez WfeC*'^'^, this reads 

k S.t. w/c^27rZ 

The control term is constructed such that the mapping in the new coordi- 
nated has j4o = as an invariant curve. In order to do this, we define the 
function $ implicitly by 

$((^) = ^ + uj + Oaw {ex'm^)), f) ■ 

The angle ^{(p) is obtained when = in Eq. (7) The expression of the 
control term is such that 

n-p) = ex'i^ + sx'iHv)) - d^w{ex'{Mv)), v). (8) 

^From the expression of w given by Eq. (4), it is straightforward to check 
that if A' is of order e then d^w = eV{K + A',(p) — eV{K, ip) is of order 
e^. Since '^{^p) — ((p + w) is of order e (again if A' is of order e) then 
ex'(^(</?)) — £x'(i^ + w) is of order . Thus, /' is of order . 
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The controlled mapping becomes 

Ao = A'o + e (x'(v^,) - x'iH^o))) + d^w {A'^ + ex'^o), V'o) - 

d^w{ex!{Hvo)).Vo). (9) 
</5o = <^o +w + aAW^(^o +^x'(¥'o)'¥'o) • (10) 

It is straightforward to sec from Eq. (9) that if = then ip^ = ^(<po) 
by definition of $ and hence Aq = . Since we assume that the mapping 
is invertible, the curve Aq = is preserved by iteration of the map. Conse- 
quently, the controlled map (5)- (6) has the invariant curve with equation 

A = ex'if). (11) 

Next, we derive an approximate control term by only keeping the order 
The expansion of /' gives the expression of /inix,2 : 

/mix,2 = ~Y^"{K) {x'iv + u:)f - e^dAV{K, (p)x'(^ + u;). (12) 

Remark 1: If we assume that w is only a function of the actions, the 
control term / is 

/(V')=x(v^^ + ^')-x($(¥'))+x'($(^))«/(c'x'($(¥')))-''^(£x'($(¥'))), (13) 

where $ is defined implicitly by ^{^p) = ip + uj + w' {ex' {^{f))). 

Remark 2: If the time step of the map is equal to r, i.e., if we consider 
controlled maps generated by 

Sc{A', if) = A'lp + tojA' + Tev{(p) + tw{A', tp) + t/(^), 

the generating function is given by 

X{^ + TUJ) - X{^) = -TV{ip). 

We define the operator 

Ht = , 

T 

and Tt as the pseudo-inverse of TLt given from "H^Fi- = Ht- The projector 
TZt is defined accordingly. Hence the solution for x is 

X={l-T,-nr)v. (14) 

The control term is 

f{v) = r-'e ix'i^ + Tco) - x'iM^))) " d^w{ex' (M^)), f), 
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where ^rif) = f + Tuj + TdAw{ex' {^t{^)) , f)- We expand the expression 
of the control term and we neglect the order t : 

/'(^) = -ex"{v')dAw{ex'{^),^) - d^w{ex\^),^) + 0{t). 

Since wx' = —v, we have 

fiip) = -w{-eTov',^)+0{T), 

where To is the pseudo-inverse of cod^. This expression of the control term 
corresponds to the one obtained by the local control of Hamiltonian flows 
in Ref.2 



3. Numerical examples 

3.1. Application to the standard map 

The standard map S is 

A' = A + e sin ip, 
(f' = (p + A' mod 27r. 

After a translation of the action A by w, the map becomes 

A' = A + e sin (p, 

ifi' = (fi + u) + A' mod 27r. 

A phase portrait of this map for e = 1.5 is given in Fig. 1. There are no 
Kolmogorov-Arnold-Moser (KAM) tori (acting as barriers in phase space) 
at this value of e (and higher). The critical value of the parameter e for 
which all KAM tori are broken is £std « 0.9716. 

The standard map is obtained from the generating function in mixed 
coordinates 

S{A', if) = A'ip + u)A + — +ecosip, 

i.e. v{lp) = cos (fi and w{A) = A^/2. The generating function x given by 
Eq. (14) is thus 

^ sin(y - uj/2) 
2sin(a;/2) ' 

The control term given by Eq. (13) is 
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Fig. 1. Phase portrait of the standard map S for e = 1.5. 



^ 




Fig. 2. Phase portrait of the controlled standard map Smix with the control term (15) 
for e = 1.5 and ui = n. The bold curve is the invariant curve created by the control term. 



where ^{(fi) — (p + ui + ex'i^if))- We notice that the equation for $ is 
invertible for s < 2sin(a;/2). 

The dominant control term is given by 



/mix,2(<p) = -£■ 



;COS^((p + uj/2) 
8sin2(w/2) ' 



(16) 



and the resulting map iSmix,2 generates the phase portrait displayed on 
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Fig. 3. 
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Fig. 3. Phase portrait of the controlled standard map 5mix,2 with the control term (16) 
for £ = 1.2. 



3.2. Application to the tokamap 

The tokamap^ has been proposed as a model map for toroidal chaotic mag- 
netic fields. It describes the motion of field lines on the poloidal section in 
the toroidal geometry. This symplectic map {A,ip) {A',(p'), where A is 
the toroidal flux and if is the poloidal angle, is generated by the function 

S{A', ip) = + H{A') - cos ^. 



It reads 



A' 

A = A' + e-^j-p^sm^, 

where q{A) = 1/H'{A) is called the q-profile. In our computation, we choose 
H'{A) = l/q{A) = 7r(2-A)(2-2A+A2)/2ande = 9/ (47r). A phase portrait 
of this map is shown in Fig. 4. 

Wc select a given value K of the action A for the localization. The 
tokamap is then generated by a function S of the form given by Eq. (1) 
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Fig. 4. Phase portrait of the tokamap for e = 9/(47r). 



with 



K 



■ cos ip, 



K+1 

w{A, if) = H{K + A)- H{K) -ojA-e 



K + A 



K 



K + A + 1 

The generating function is given by Eq. (14) : 

K sin(</j — w/2) 
" 2{K + 1) sin(w/2) ■ 

The expression of /mix,2 is given by Eq. (12) : 



K+l 



cos if. 



mix, 2 



e^K cos((^ + tj/2) 
2(if+l)2 sin(tj/2) 



cos ip 
K + l 



KH"{K) 



cos((^ + w/2) 



sin(a;/2) 



For K = 1/2 the controlled tokamap is 



A = A' + e 



77 Sin cp — — 

l + A' ^ 9 

1 e 



2sm{2ip + a) Il7r2(^ + 2a 



sma 



+ 



64 



sm Q 



mlCos</j, 



(17) 
(18) 



q{A') {1 + A'f 

where a = ll7r/32 We notice that this control term is the same as the 
one obtained by performing the global control and then by expanding the 
control term around a given value of the action. The main advantage here 
is that the whole scries of the control term can be computed which was not 
the case with the global control. 
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Fig. 5. Phase portrait of the controlled tokamap (17)-(18) for e = 9/(47r). 
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